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ABSTRACT 

A  method  of  determining  acoustic  power  output  of  undent 
sound  sources  is  presented  employing  reverberation  techniques 
similar  to  those  employed  in  air. 

Two  possible  approaches  were  taken „  One  following  the  id* 
of  R0  W„  Young,  which  includes  the  effect  of  the  direct  wave  and 
by  use  of  a  correction  factor  eliminates  its  effect  on  the  result- 
ant power,  which  requires  knowledge  of  a  steady  state  average 
pressure,  a  decay  rate  and  an  average  absorption  coefficient „ 

The  second  follows  the  idea  presented  in  Fundamentals  of 

11 

Acoustics  by  Kinsler  and  Frey  which  requires  knowledge  of  only 

a  steady  state  average  pressure  and  decay  rate0 

Both  methods  yield  results  consistent  to  within  one  decibel 
of  a  standard  acoustic  power  as  determined  by  pulse  techniques „ 

A  test  setup  is  proposed  for  use  of  this  technique  with 
comments  given  as  to  desirable  characteristics  of  equipment « 
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SYMBOL  TABLE 

a  Sabine  absorption  in  Sabins,?  equals  Sa 

a  average  Sabine  coefficient  for  a  room 

a.  random  incidence  energy  absorption  coefficient  of  a 

particular  surface  of  area  S. 

l 

°<  average  random  incidence  energy  absorption  coefficient 

c  speed  of  sound  in  the  medium  in  the  room  or  enclosure 

D  decay  rate  of  sound  pressure  level  db/sec 

£  average  acoustic  energy  density  equals  P  jpo. 

E  RMS  voltage 

f  center  frequency  of  noise  band 

I  acoustic  intensity 

K  constant 

P  acoustic  pressure(rms ) 

£>  density  of  the  medium 

R  equivalent  parallel  resistance  of  hydrophone 

S  total  surface  area  of  an  enclosure 

S.  area  of  a  specific  side  of  the  enclosure 
1 

T  reverberation  time  in  seconds  equals  bO/D 

V  volume  of  the  enclosure 

W  acoustic  power  output  of  the  source  in  watts 
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1.   Introduction. 

The  determination  of  absorption  coefficients  of  acoustical  materials 
has  attracted  considerable  interest  and  effort  over  the  years. 

Recently  R.  W.  Young  published  a  review  article  in  this  field  which 
proposed  modification  to  present  engineering  practices „   The  paper  also 
presented  a  method  for  determining  the  sound  power  output  of  a  transducer 
by  measurement  of  decay  rate,  a  steady  state  sound  pressure  and  an 

average  absorption  coefficient  in  the  enclosure  under  consideration. 

2 
R.  W.  Case  and  E.  W.  Vahlkamp  made  measurements  of  power  through 

decay  rates  as  proposed  by  Young  but  were  unable  to  achieve  consistent 

agreement  with  power  measured  by  pulse  techniques  „ 

3 
Kinsler  and  Frey  also  suggest  a  second  or  alternative  equation 

for  the  calculation  of  acoustic  power  from  measurement  of  decay  rate 

and  a  steady  state  pressure. 

The  purpose  of  this  investigation  is  to  test  the  validity  of 

these  theories  on  power  measurements  as  applied  to  a  reverberation 

tank  and  a  water  medium. 


R.  W.  Young,  Sabine  Reverberation  Equation  and  Sound  Power 
Calculations,  J.A.S.A.,  31,  pp  912-921  July  1959 

p 
R.  W.  Case  and  E.  W.  Vahlkamp,  Investigation  and  Application  of 

Reverberation  Measurements  in  Water,  USNPGS,  Thesis,  I962 

L.  E.  Kinsler  and  A.  R.  Frey,  Fundamentals  of  Acoustics 
Second  Edition,  John  Wiley  &  Sons,  I962 


2.   Theory. 

All  the  contributions  to  architectural  acoustics  betore  the  work 
of  W.  C 0  Sabine  were  of  a  qualitative  nature „   Sabine  was  the  first  to 
commence  a  comprehensive  and  quantitative  study  of  the  subject  and^ 
since  science  and  engineering  are  basically  quantitative ^  it  may  be 

fairly  stated  that  the  science  of  architectural  acoustics  had  its 

k 
beginning  with  Sabine  „   In  his  paper  in  1900  ,,  Sabine  set  forth  in 

simple, but  comprehensive, language  the  requirements  for  good  hearing  in 

any  auditorium  as  follows: 

It  is  necessary  that  the  sound  should  be  sufficiently  loud; 

that  the  simultaneous  components  of  a  complex  sound  should 

maintain  their  proper  relative  intensities;  and  that  the 

successive  sounds  in  rapidly  moving  articulation,,  either 

speech  or  music,  should  be  clear  and  distinct^,  free  from 

each  other  and  from  extraneous  noises „   These  three  are  the 

necessary,  as  they  are  the  entirely  sufficient^,  conditions 

for  good  hearing o 

It  seemed  apparent  to  Sabine  that  the  third  of  these  factors  was 

the  most  important  one  in  affecting  the  acoustical  quality  of  a  room 

and  he  therefore  devoted  a  good  number  of  years  to  a  quantitative 

study  of  the  growth  and  decay  of  sound  in  an  enclosure.   From  this 

study  he  formulated  the  following  empirical  equation^ 

T  -  H" 


W„  Co  Sabine,  Collected  Papers  on  Acoustics^,  Harvard 
University  Press,  1922 


where:  T  = 

k  = 

V  = 

S  = 

a  = 


reverberation  time 

constant 

volume  of  the  enclosure 

total  surface  area  of  the  enclosure 

average  Sabine  absorption  coefficient  of 

the  enclosure 
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Several  authors  '  '   set  forth  a  derivation  following  the  method 

o 
developed  by  Jaeger  based  on  the  following  assumptions: 

1.   A  random  or  diffuse  distribution  of  the  flow  of  sound  in 

the  room  is  required  and  2.  a  continuous  absorption  of  sound 

by  the  boundaries  of  the  room  occurs.   These  conditions  are 

quite  completely  fulfilled  in  highly  reverberant  rooms.   It 

is  to  be  expected  therefore  that  these  formulas  would  be 

satisfactory  for  the  practical  calculations  of  reverberation 

in  live  rooms  but  that  they  would  lead  to  only  approximate 

results  in  dead  rooms.   To  say  a  room  is   live  we  mean  that 

the  absorptive  power  is  small,  the  rate  of  decay  is  slow  and 

hence  the  reverberation  time  is  long. 

The  assumptions  above  give  rise  to  the  following  differential 

equation  with  regard  to  the  conservation  of  energy: 


Rate  of  increase 
of  energy  in  room 


Rate  of  emission 
of  energy  from 
source 


Rate  of  absorption 
of  energy  by  walls 


(2) 


V.  0.  Knudsen,  Architectural  Acoustics,  John  Wiley  and  Sons,  Inc., 
1932 

P.  M.  Morse  and  R.  H.  Bolt,  Sound  Waves  in  Rooms,  Rev.  Mod.  Phys., 
16,  p  76,  Jan  ISkk 

L.  E.  Kinsler  and  A.  R.  Frey,  op.  cit. 

A.  Jaeger,  Zur  Theorie  des  Nachalls,  Akad,  Wiss,  Wein  2A, 
120,  p  613,  1911 


Let  £.  be  the  average  acoustic  energy  density^  V  the  volume  of  the 
enclosure,  and  W  the  power  output  of  the  source.   Then  the  rate  of 
absorption  of  energy  by  the  walls  can  be  shown  to  be 

£^A  (3) 

where  c  is  the  speed  of  sound  and  a  is  the  total  absorptive  power  of 
a  live  room  which  is  found  from 

a  =•  2  *i^t  (k) 

where  a.  are  the  respective  absorption  coefficients  representing 
the  fraction  of  randomly  incident  energy  absorbed  by  the  different 
materials  of  areas  S.  forming  the  interior  walls  of  the  room  as 
well  as  any  other  absorbing  surfaces. 

The  rate  at  which  the  energy  increases  in  the  room  is  given  by 

V^  (5) 

Substitution  of  equations (3)  and  (5)  into  (2)  yields 

VH  * w "  H*  (6) 

Assuming  the  sound  source  to  have  been  started  at  t  =  0  the  solution 
of  equation  (6)  is 


-  act, 


>  I  -     Oily 


This  can  be  written  in  terms  of  the  mean  square  acoustic  pressure 
by  use  of  the  relation  that 


t    -  ^c?  (8 


This  results  in 


.9) 


If  the  source  is  left  on  until  steady  state  is  reached^,  ice„  t  s  OO 
equation  (9)  reduces  to 

r  s  — "f"*  (10  j 

or 

The  differential  equation  governing  the  decay  of  uniformly 
distributed  diffuse  sound  in  a  live  room  is  obtained  by  letting 
W  =  0  in  equation  (6).   If  we  assume  the  source  is  turned  off  at 
t  =  0  and  letting   £0  represent  the  assumed  uniformly  distributed 
energy  density  at  this  instant  then  as  t  increases 

£  -  £c  e""^7-  (12) 


or  in  terms  of  intensity 

-3ct 


Applying  the  operator  lOlog  to  both  sides  of  the  equation  it  becomes 

act 
A  XL   «  »o  1°^ 


,   ©  "W    -  '.o*7  act 


where  AIL  =  lOlog  I/I  represents  the  change  in  intensity  level  in 
decibels.   It  can  now  be  said  that  the  intensity  level  in  a  live 
room  and  correspondingly  the  sound  pressure  level  decreases  with 
elapsed  time  at  a  constant  decay  rate  D  in  decibels  per  second  given 

D  -    g (1*0 

In  accordance  with  the  original  idea  of  Sabinej  let  us  define 
reverberation  time  T  as  the  time  required  tor  the  level  of  the  sound 


in  the  room  to  decay  by  60  db .   Then 

T  -  D   ~   a  c  (15} 

which  is  the  form  of  Sabine  s  emperical  equation „ 

It  can  be  seen  however  that  under  dead  room  conditions^  i„et 
a  very  large,  approaching  unity,  the  predicted  reverberation  time  by 
equation  (15 )  is 

-r     ^  2  V 

T  =    Sc  (16) 

which  is  not  true  since  we  know  that  for  a  perfect  absorber  the 
reverberation  time  is  zero. 

Eyring  produced  the  first  important  modification  to  Sabine's 
equation.   His  derivation  was  based  on  image  sources  <,  all  of  which 
came  into  existence  when  the  real  source  starts.   The  growth  of  the 
intensity  is  then  an  accumulation  of  successive  increments  from  the 
true  source,  from  the  first  order  or  single  reflection  images^  from 
second  order  or  double  reflection  images  and  so  on  until  all  sources 
of  any  appreciable  strength  have  made  their  contribution   At  this 
time  steady  state  is  reached  and  the  rate  of  absorption  and  rate  of 
emission  of  sound  are  equal,,   Similarly  when  the  true  source  is 
stopped,  all  the  images  are  stopped.   The  decay  in  energy  density  in 
the  room  therefore  results  from  the  successive  losses  of  acoustic 
radiation^  first  from  the  true  source,  then  from  the  first  image  and 
so  on  until  all  the  images  have  radiated  their  energy  into  the  room. 
From  this  we  see  that  the  duration  of  audibility  is  equal  to  the  time 
required  for  sound  to  arrive  at  the  room  from  the  order  of  images  which 


Q  HI      HI 

C.  F „  Eyring,  Reverberation  Time  in  Dead  Rooms 
J.A.S.A.,  1,  pp  217-241,  Jan.  1930 
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is  so  far  removed  from  the  room  that  all.  images  beyond  this  order  will 
contribute  enough  energy  to  be  60  db  down  from  the  steady  state  energy, 
It  can  be  seen  that  the  Sabine  assumption  of  uniform  distribution 
of  energy  and  random  flow  of  energy  are  fulfilled  but  the  decay  is  not 
continuous  but  occurs  in  steps  resulting  in  a  greater  absorption 
during  the  same  interval  of  time  and  a  more  rapid  decay  rate 0  This  is 
the  essential  difference  in  Eyring  s  equation  and  when  used  in  dead 

rooms  gives  results  superior  to  the  Sabine  equation,, 

10  s 

Norris   suggests  a  method  for  deducing  Eyring  s  equation  with 

the  result  being  the  following  equation; 

or  in  terms  of  decay  rate 

^    -  |.0?7  cS  \w  (\  -«0 

■J)  *   _ (18) 

where  5  =  (£  c»i.Si  V^  $L  (19) 

ii 
It  should  be  noted  the  Eyring  equation  uses  an  average 

absorption  coefficient  .   It  is  assumed  that  the  walls  are  uniformly 

covered  with  material  of  uniform  absorption  or  that  the  material  is 

sufficiently  well  spread  over  the  walls  so  that  an  average  value  can 

be  taken. 

When  o<  is  small  it  can  be  seen  that  Eyring  s  equation  reduces 

to  Sabine's  equation  which  is  to  be  expected  for  a  live  room. 

However,  it  should  be  noted  that  the  Eyring  form  of  the  equation 

is  seriously  in  error  for  very  nonuniformly  placed  absorbents.  An 


10, 


R.  F.  Norris,  A  Derivation  of  the  Reverberation  Formula, 
Appendix  II  in  Architectural  Acoustics  by  V.  0.  Knudsen, 
John  Wiley  and  Sons,  1932. 


example  of  this  would  be  an  enclosure  in  which  one  wall  is  highly 

absorbent  and  the  others  are  highly  reflective  „ 

11  12 
Millington  and  Sette   '    proposed  another  form  of  averaging  the 

absorption  coefficient.   Their  equation  is: 

The  main  point  of  difference  is  this:   Norris  -  Eyring  s  theory 
assumes  that  the  energy  in  the  room  resumes  uniform  distribution  after 
each  set  of  incidences  during  the  discontinuous  decay  processes^ 
Millington  and  Sette  follow  the  course  of  a  bundle  of  rays  through 
many  reflections  and  assume  that,  on  the  average,,  a  particular  ray 
will  strike  a  given  surface  a  number  of  times  proportional  to  its 
area.   Both  forms  assume  the  Sabine  geometric  conditions  but  the 
averaging  is  obtained  differently,,   Millington  and  Sette  take  a  geo- 
metric mean  whereas  Norris  -  Eyring'  s  is  an  arithmetic  mean,,   The 
serious  defect  of  the  Millington  and  Sette  form  is  that  it  predicts 
T  =  0  if  any  surface,  no  matter  how  small  in  area^  is  a  perfect 

absorber . 

13  i 

Fitzroy   proposes  another  modification  to  Eyring  s  equation 

which  yields  reasonable  results  for  rooms  in  which  the  absorption 

coefficients  differ  widely.  This  formula  is  different  as  follows: 


UG.  Millington,  J.A.S.A.  k,    pp  69-82,  1932 

12W.  J.  Sette,  J.A.S.A.  k,    pp  193-210,  1933 

13D.  Fitzroy,  J.A.S.A.  31,  pp  893-897,  1959 
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in  the  case  of  a  rectangular  room,  three  different  calculations  are 
made  by  means  of  the  Eyring  equation  but  the  average  absorption  is 
changed  for  each  pair  of  boundaries „   Ratios  are  established  relating 
each  pair  of  boundary  areas  to  the  total  room  area0 
This  equation  can  be  written  thus  % 


T  = 


SL-S  w(i-*v>P  S1-S  W(v--ay\r  ![-$>  wv(»-aij  (21) 


Where;   x  =  total  side  wall  area 

y  =  total  ceiling  and  floor  area 

z  =  total  end  wall  area 

S  =  x+y+z 

ax  =  average  absorptivity  in  x  area 

a  =  average  absorptivity  in  y  area 

a  =  average  absorptivity  in  z  area 


It  should  be  noted  here  that  care  must  be  taken  to  determine 

whether  the  room  is  live  or  dead  and  to  determine  which  equation  will 

Ik 
best  fit  the  given  geometry.   To  quote  Knudsen   , 

The  approximate  theories,  when  used  with  caution  and  under- 

standing,  have  served  satisfactorily  for  practical  purposes 

of  acoustical  designing  and  they  will  continue  to  do  so 

until  they  are  superseded  by  more  exact  theories „ 

The  application  of  room  acoustics  to  a  water  medium  seems  to  be 

the  solution  to  one  of  the  problems  that  has  arisen  in  underwater 

measurements . 


V.  0,  Knudson,  Recent  Developments  in  Architectural 


Acoustics,  Rev.  Mod.  Phys,  6,  p  3j>  Jan  193^ 


Test  tanks,  which  have  simulated  free  field  conditions  for 
small,  high  frequency  transducers  of  the  past,,  are  no  longer  good 
free  field  approximations  for  the  large,,  low  frequency  transducers 
of  today.   Other  problems  have  also  arisen „   The  solution  or  partial 
solution  to  these  problems  seems  to  lie  in  the  use  of  a  reverberation 

chamber . 

15 
In  Kinsler  and  Frey   an  equation  is  given  for  measuring  acoustic 

power  output  of  sound  sources  in  a  reverberant  chamber,, 

If  we  take  equation  (11)  and  substitute  for  a  from  equation  (14)  the 

following  results: 

T*VD 


W  = 


22 


4  I  l.otfi},*  c* 

It  is  interesting  to  note  that  this  equation  is  independent  of 
a  and  does  not  require  knowledge  of  the  absorption  coefficients  of 
the  room  as  long  as  it  is  live.   If  we  rewrite  our  differential 
equation  (6) 


18. 


m    * 


ac 


=  W 


(23; 


and  substitute  a  =  -Sln(l~*)  after  Eyring^  we  find  the  solution 
to  be 


4W 


^  =  :^tW'-*H 


4  V 


2k 


or  in  terms  of  average  pressure 


,25 


Kinsler  and  Frey,  op.cit.,  pp  ^36 

10 


at  steady  state  or  t  =  oO  f   we  find 


substituting  equation  (l8)  for  decay  rate  we  find 

4(1.027)^0.^ 
Which  is  identical  to  equation  (22)  above. 

Young   presents  several  different  forms  of  power  equations  to 
obtain  power  output  of  a  transducer  by  measurement  of  decay  rate,, 
steady  state  sound  pressure  and  an  average  absorption  coef f icient „ 
One  of  the  forms  presented  suitable  for  use  in  a  reverberant 
chamber  is 


28 


This  equation  is  derived  from 


W  =  ^S    *  (29) 

4f>C 

which  has  been  arrived  at  not  by  ignoring  the  direct  sounds,  but  by 
adding  it  to  the  reflected  sound.   By  substitution  of 

5    «    i  -  e  (30) 

and 

S     a  V"^  -  (31) 

equation  (28)  above  results. 

It  should  be  noted  that  for  small   d   this  expression  reduces  to 
equation  (22). 

R.  W.  Young,  op.  cit.,  pp  916-918. 

11 


3.   Description  of  tank  and  lining  material , 

A  small  free  standing  tank  (l„82  x  1.22  x  variable  depth  to  a 
maximum  of  0,75  meters)  was  used  in  making  the  reverberation  measure- 
ments.  It  was  constructed  from  0oQQk   meter  thick  black  iron  sheet 
and  set  upon  0.02  meter  thick  rubber  pads  resting  on  0d5  meter  wooden 
blocks . 

An  unlined  tank,  its  interior  painted  with  a  glossy  enamel^,  was 
first  investigated.   The  second  condition  under  which  the  reverberation 

phenomenon  was  studied  was  with  a  few  blocks  of  cone  lattice  structure 

17 
aluminum  loaded  butyl  rubber  lining   „   These  blocks  are  0.2  x  0.^-1 

meter  and  were  used  in  varying  numbers  from  one  to  twelve  which  re- 
sulted in  a  variable  decay  rate. 

It  should  be  noted  that  the  maximum  decay  observed  with  12  blocks 
at  50Kc  results  in  an  a  of  O.I65. 


17 

A.  Heller,  Underwater  Anechoic  Tank  Linings,  Navord  Report  29^9 

(Confidential),  U.  S.  Naval  Ordnance  Laboratory^   White  Oak, 

Maryland,  pp  19-22,  10  Nov  1953 
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k .   Description  of  equipment  and  technique. 

Figure  1  shows  the  block  diagram  of  the  equipment  used  in  rever- 
beration measurements.,   The  low  level  noise  signal  was  passed  through 
a  bandpass  filter  having  the  following  characteristics^  lower  3^b 
point  at  „78f  ,    upper  3db  point  at  1.28f  and  l8db  per  octave  slope „ 
The  General  Radio  1391  Pulse  and  Time  Delay  Generator  was  triggered  by 
the  square  wave  generator  which  controlled  the  repetition  rate  of  the 
noise  pulses .   When  repetition  rates  lower  than  one  per  second  were 
required,  a  battery  and  hand  key  were  substituted  for  the  square  wave 
generator.   In  this  way  the  noise  signal  could  be  turned  on  and  off 
at  a  rate  adjustable  to  allow  full  buildup  and  decay  of  the  sound 
field.   The  filtered  noise  signal  was  then  amplified  to  a  value  of 
3.3  volts  RMS  which  was  maintained  throughout  the  investigation  and 
transmitted  into  the  water  by  an  LC-32  hydrophone ,   The  receiving  portion 
of  the  setup  consisted  of  a  similar  LC=32  hydrophone  followed  by  a  linear 
amplifier,  logarithmic  amplifier  and  Memoscope.   For  low  decay  rates^, 
up  to  a  maximum  of  approximately  200  db/sec,  the  Bruel  and  Kjaer  recorder 
was  useful. 

The  Memoscope  vertical  scale  was  calibrated  in  db  to  facilitate 
reading  the  decay  rate  in  db/sec .   The  storage  facility  of  the  Memoscope 
permitted  retracing  the  pattern  and  then  by  use  of  a  Polaroid  camera  a 
permanent  record  was  made . 

The  same  equipment  was  used  to  make  steady  state  pressure  measure- 
ments employing  a  continuous  signal  and  the  Fluke  Model  910A  true 
RMS  voltmeter. 
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The  equipment  for  pulse  technique  measurements  is  essentially  the 
same  as  Figure  1  with  the  substitution  of  a  Western  Electric  17B 
oscillator  in  place  of  the  noise  source  and  filter  „   Measurement  of 
input  and  output  voltages  were  made  by  suitably  calibrating  the 
Memoscope  with  a  known  voltage., 

The  filtered  noise  signal  resulted  in  a  considerably  diminished 
standing  wave  ratio  of  approximately  2db  as  compared  to  approximately 
30db  for  a  sine  wave. 

It  should  be  noted  that  the  tank  linings  were  not  uniform  but 
the  a  is  such  that  the  tank  is  live,,  i.e.  a<  0.2  and  Sabine's 
equation  and  theory  applies.   The  actual  value  of  a  is  not  important 
to  the  calculations  since  decay  rates  are  used. 

The  testing  procedure  consisted  of  taking  reverberation  measure- 
ments with  no  lining,  six  blocks  lining  and  twelve  blocks  lining  with 
the  addition  of  a  balloon  as  noted  below.  These  conditions  were 
repeated  over  a  period  of  several  days  which  required  removal  of  the 
blocks  from  the  tank  in  order  to  start  a  new  series  of  measurements. 
The  initial  placement  of  the  blocks  was  in  a  random  manner.   Success- 
ive tests  were  made  with  the  blocks  placed  in  the  same  approximate 
position  so  as  to  minimize  the  effect  on  decay  rate  by  variations  in 
placement  of  the  sound  absorbing  material.   It  is  thought  that  this 
step  is  really  not  necessary  since  a  is  <C  0.2  but  this  thought  was  not 
fully  investigated. 

Experiments  were  made  with  various  diffusers  to  eliminate  or 
reduce  the  direct  wave.   The  end  result  was  that  a  balloon  filled 
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with  air  and  weighted  to  hold  it  in  place  served  best.   When  placed 
in  line  between  the  transducers  the  amplitude  of  the  direct  wave 
was  reduced  one  order  of  magnitude.   The  balloon  adds  absorption 
and  the  result  is  an  increase  in  decay  rate  from  97  db/sec„  in  the 
unlined  tank  without  the  balloon  to  192  db/sec0  in  the  unlined  tank 
with  balloon. 

As  more  absorptive  material  is  added  the  reverberant  field 
becomes  less  intense  and  less  diffuse,  resulting  in  a  larger  standing 
wave  ratio.   The  balloon  also  serves  to  break  up  this  standing  wave 0 
For  this  reason  the  balloon  was  left  in  the  tank  when  investigating 
Young  s  equation,  but  was  moved  out  of  the  path  of  the  direct  wave  „ 

The  balloon  could  be  described  as  being  the  shape  of  an  inverted 

pear  with  a  length  of  25  cm  and  a  diameter  of  12  cm. 

18 
According  to  Allred  and  Newhouse   there  is  a  slight  disagreement 

between  the  mean  free  path  obtained  by  use  of  the  equation  kv/S 

and  the  value  obtained  by  Monte  Carlo  methods „   In  the  situation  under 

study  it  is  noted  that  our  tank  1.8X1 .2X0.5  can  be  scaled  to 

10X6.66X2.78.   By  comparing  this  with  Allred  and  Newhouses  10X6X2  we 

see  the  correction  factor  is  1.032  or  a  3»2$  error  is  made.   This  is 

considered  negligible  in  the  case  under  study  but  for  work  of  the 

most  precise  nature  should  be  considered. 


1  Q 

Allred  and  Newhouse,  Applications  of  the  Monte  Carlo 
Method  to  Architectural  Acoustics,  J.A0S0A0,  30^  pp . 
903-90^,  Oct.,  1958. 
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5c   Results. 

From  the  theory  presented  previously  we  must  formulate  wor 
equations  to  take  into  account  the  directivity  and  receiving  response 
of  the  LC-32  hydrophones  used. 

The  directivity  ratio  of  a  hydrophone  is  given  by 

t>  =  ^=1   -  ^L.  (32) 

where;   D  =  directivity  ratio  subscripted  with  a  t  or  r  to 

indicate  whether  transmitting  or  receiving  direct- 
ivity. 

I   ,  =  the  intensity  along  the  reference  axis  at  some 
re  f 

distance  r  from  the  source. 

I    =  the  average  intensity  defined  by  I    s  W/4TTr~ 
ave  J  J      ave    ' 

where  W  is  the  total  acoustic  power  radiated 

by  the  source. 

The  same  notation  applies  to  P,,  the  RMS  pressure  derived  from 

2 
the  relation  I  =  P  /p  C     . 

By  pulse  technique  it  is  possible  to  measure  a  standard  acoustic 
power  output  and  from  this,  knowing  the  input  voltage  and  equivalent 
parallel  resistance  of  the  hydrophone^  the  efficiency  can  be  calcu- 
lated.  The  resultant  power  output  and  efficiencies  are  then  compared 
with  those  observed  by  reverberation  technique. 

Pulse  technique  power  is  given  by 

W  -       +^lV-"  (33) 

In  the  method  used  a  voltage  proportional  to  the  axial  pressure  is 
measured.   Hence  the  average  pressure  is  found  by 


e 


*       £*  (34) 


AVE        p 
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where  E  is  the  RMS  voltage  measured 

M  is  the  hydrophone  receiving  response  in  volts  per 
newtons/square  meter  =  E/P 

Thus  our  working  equation  becomes 

W  =   7/'f*  (35, 

where  D   is  the  transmitting  directivity  of  the  LC~j2  hydro- 
phone 

In  the  reverberation  technique  the  power  is  given  by  equation 

(22) 

^  -   4(i.o87Vci  (36, 

In  this  technique  the  voltage  measured  is  proportional  to  the 
average  pressure  but  since  the  response  of  the  LC-32  hydrophone  is 
given  with  respect  to  the  axial  pressure,,  the  tollowing  expression 
results 

A  discussion  or  directivity  ratios  is  given  in  Appendix  IQ 
Sample  Calculations; 

A,   Pulse  Power  at  JOKc  using  equation  (35)* 


f>c  n1^ 


r  =    „31fc>  meter 
.3 


p  =  998  Kg/r 
c   =   1485.8  m/s   (at  2l.5°C; 
^dk=  -105.08   db  l/M2  =   3.215   x   108  Newto^/meter1*  volt2 

Dt-1.63 

E   =   l.W  x    10   J   volts   RMS 
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^   .  34.  *f  milliwatts 


B„   Reverberation  Power  at  30  Kc  using  equation  (37] 
no  lining,  no  balloon. 


V\x  4Ci.og-7)^c" 


V  =  1.08  meter3 
D  =1.60 
D  =  97  db/sec 
E  a  8A5  x  lO"3  volts  RMS 


.    Q.oyyt-4oM3.Zl*,*/o*)(<H)(g,*£'y<o-3) 
<l-(l.«>*'?)(W)(l*»-ttr.&)* 

=   >o3  milliwatts 


C.   Reverberation  Power  at  30  Kc  using  equation  (28) 
modified,  12  blocks,  balloon  not  in  path  of  direct 
wave . 


D  =  1675  db/sec. 

S  =  7.32  meter^ 

E  =  2.14  x  10" -*  volts  RMS 

from  equation  (lU) 


a  =■ 


I?  V 


.ofrT  Sc 


1.0*7  (7.  ii)(ilf8i'  i) 


I'd 


w  ^ 


E'DVD, 


t"V  4(i.oSf7)/3Cx    |_         5 


NH 


(2./4x/o"*)  (iUns)Ll^i)(i-J>o)ll^i<:rtfo2)\  ■   -  e 


*M  I .  »&7)  (<?<?&)  (i4*6-.SV 


4t±     milliwatts 
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Comparison  of  these,  i.e.  Reverberation  power  against  Pulse 

power  is  not  readily  obvious  due  to  the  difference  in  true  input 

power  for  the  same  voltage  input  0   It  is  not  to  be  expected  that 

loading  of  the  medium  will  be  the  same  since  in  the  case  of  the  pulse 

technique  the  transducer  sees  essentially  a  free  fields  but  sees  a 

standing  wave  pattern  under  reverberation  conditions „ 

Measurements  of  R  were  made  under  these  conditions  with  the 
P 

results  being  shown  in  Table  II-2,  Appendix  II c 

We  can  best  compare  the  power  by  normalizing  the  pulse  power 

to  an  equivalent  reverberation  power »  The  procedure  is  as  follows t 

Let  W  =  the  output  acoustic  power  measured  under  pulse  technique 
P 

with  an  input  voltage  E.  volts,  transducer  equivalent  parallel 
resistance  R   and  efficiency  \\      .   Then  we  can  say 
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Similarly  for  the  reverberation  case  let  W  s  the  measured  acoustic 

R 

power>  E2  the  input  voltage  (held  constant^  E^  »  E  )9    R   the  parallel 
resistance  and  v^^the  efficiency.   Then 


'R  - 


*f*         ftT3t 


Since  we  are  considering  the  same  transducer  at  the  same  center 
frequency  v\t    =  v^  .   Thus 

In  our  case  at  30  Kc  we  find 

R  .  =  6.55  K  ohms 
Pi 

R  _  =  5.76  K  ohms 
p2 

Thus  we  can  say  that  the  power  measured  by  pulse  technique  being 

0.364  milliwatts  should  be  normalized  to 

,2,44  *  J^A —   _-   .  4.1  *f  milliwatts 

Let  us  call  the  new  value  the  standard  or  true  value  of  the  power. 

Now  we  can  compare  our  reverberation  power  to  the  standard 
value.   Table  1  is  such  a  comparison. 

Figures  2  thru  8  show  typical  decay  curves  for  30KC  under  the 
various  conditions  studied. 

A  comment  should  be  made  about  the  notation  used  in  Table  1  and 
in  the  figures.   Balloon  in  path  means  the  balloon  was  placed 
between  the  transducers  blocking  the  path  of  the  direct  wave  and 
reducing  its  amplitude  by  approximately  one  order  of  magnitude. 
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.     .  51  I/? 

Equation  (37)  is  used  for  these  calculations.    Balloon  not  in  path 
means  the  balloon  was  in  the  tank  but  not  between  the  transducer 
the  path  of  the  direct  wave.   Equation  (28 )s   modified  as  shown  in 
sample  calculation  C,  is  used  for  these  calculations. 
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Table  1 

Acoustic  Power  Output  (milliwatts)  for  3  3V  RMS  Input 

Center 

Frequency      10     15     20    25     30         4Q    45 

Kilocycles 

Pulse  Tech. 

Measured       .0102   ,0381   .105   .1775   .364   .59I   1,01   1.88  2 .98 

Pulse  Tech. 

Standard1      .0123   .0474   .131   .221    .414   .669   1.14  2.02   3.12 

No  lining 

No  balloon     .0127   .0458   .117   .195    M2      „684   1.05   2.05   3.00 

No  lining, 

balloon  in     .0138   .0495   .115   .194   ^39  «763  1.18  2.37  3.32 

path 

No  lining, 

balloon  not    .0128  .0558  .149   .256  .512  .835  I.28  2.38  3.61 

in  path 

6  blocks, 

balloon  in     .0118   .0417   .121    .199   .394   .615   1.05  1.90  2.88 

path 


6  blocks, 

balloon  not     .0110   .0446   .131    .211   .476   .800  1.21  2.20 

in  path 

12  blocks, 

balloon  in     .0106   .0451   .127    .198   .408   .642   1.04  2.09  2.89 

path 

12  blocks, 

balloon  not     .0129   .0481   .110    .194   .425   .667   1.04  2.02   -   j 

in  path 

This  is  pulse  power  when  normalized   for  differences  in  R  as 
noted  in  text. 
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Vertical  Scale: 
l)    db/division 

Horizontal  Scale  : 
l)0   msec/division 

E  =  (-  . 

D  =  97  db/sec 

W  =  .UlU  milliwa 


Figure  2 
Photograph  of  Decay  Rate;  No  lining  No  Balloon,  30Kc 


Vertical  Scale; 
5  db/division 

Horizontal  Scale  : 
20  msec/division 

E  =  6.2  millivolts  RMS 
D  =  192.5  db/sec 
/J  =  .4^0  mi  1  liwat  ts 


Photograph 


Figure  3 

Path, 


Vertical  Scale: 
5  db/division 

Horizontal  Scale  : 
20  msec/division 

E  =  '  .     illivolts  RMS 

D  =  192  db/sec 

W  =  .515  milliwatts 


1  re  k 
Photograph  of  Decay  Rate;  Ao    Lining,  Balloon  not  in  Path,  jOKc 


Vertical  Scale: 
5  db/division 

Horizontal  Scale  : 
5  msec/division 

E  =  2.6b  millivolts  RMS 

D  =  9*4-0  db/sec 

W  =  o9^  milliwatts 


Figure  b 

Photograph  of  Decay  Rate;  Six  Blocks  of  Lining,  With  Balloon  in  Path, 
jOKc 


Vertical  Scale: 
5  db/division 

Horizontal  Scale: 
5  msec/division 

E  =  2.^5  millivolts  RMS 

D  =  ^50  db/sec 

W  =  .472  milliwatts 


Figure  6 

Photograph  of  Decay  Rate;  Six  Blocks  of  Lining,  Balloon  not  in 
Path,  30Kc 


Vertical  Scale; 
5  db/division 

Horizontal  Scale: 
2  msec/division 

E  =  1.98  millivolts  RMS 

D  =  1750  db/sec 

W  =  .U06  milliwatts 


Figure  7 

Photograph  of  Decay  Rate;  12  Blocks  of  Lining,  With  Balloon  in  Path, 
jOKc 


Vertical  Scale: 
5  db/division 

Horizontal  Scale  : 
2  msec/division 

E  =  2.15  millivolts  RMS 

D  =  1650  db/sec 

W  =  .k2^   milliwatts 


Figure  8 

Photograph  of  Decay  Rate;  12  Blocks  of  Lining,  Balloon  not  in  Path 
JOKc 


6.   Conclusions. 

From  the  results  it  can  be  seen  that  under  conditions  of  equal 
efficiency  the  powers  measured  are  well  within  one  decibel.   Thus  the 
reverberation  technique  could  be  used  to  measure  acoustic  power 
output  with  reasonable  accuracy . 

The  reverberation  technique  appears  to  have  certain  advantages 
over  the  pulse  technique  or  other  techniques.   These  ares 

1.  Relatively  simple  and  quick. 

2.  The  output  voltage  can  be  measured  directly  by  use  of  a 
true  rms  meter  and  does  not  require  complicated  gating  setups  as 
required  in  some  pulse  setups. 

3.  Gives  results  within  1  db  of  pulse  powers  which  in  them- 
selves have  about  1  db  accuracy. 

The  disadvantages  of  such  a  system  are; 

1.  The  slope  of  the  decay  curve  must  be  evaluated. 

2.  The  source  of  the  noise  power  must  be  controllable 9    i.e. 
it  must  be  capable  of  being  shut  off  when  desired  „   If  the  source 
is  not  controllable,  a  secondary  source  with  the  same  frequency 
can  be  used  to  determine  the  decay  rate  for  the  enclosure „   The 
source  under  test  can  then  be  used  to  determine  the  average  steady 
state  pressure. 

One  possible  situation  that  arises  is  in  the  measurement  of 
machinery  noise  output  of  a  submarine  or  surface  ship  in  an  enclo= 
sure  such  as  a  drydock.   By  measuring  an  average  pressure^,  either 
by  rotating  the  hydrophone  or  by  measurements  in  varied  locations^ 
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and  measuring  the  reverberation  time  equation  (22)  can  be  used  to 
calculate  the  power.   This,  of  course,,  includes  the  direct  wave* 
Equation  (28)  can  be  used  if  an  average  absorption  coefficient  is 
calculated o   It  should  be  noted  that  the  difference  between  the  power 
calculated  by  equation  (28)  is  different  from  that  calculated  by 
equation  (22)  by  the  factor  (1-2   )  /a  which  for  a  0„2  results  in  less 
than  IQfjo   error.   Thus  particular  care  need  not  be  taken  to  ensure 
elimination  of  the  direct  wave,  provided  the  lest  tank  is  live 

This  technique  could  also  be  used  in  determining  the  effect  of 
noise  reduction  treatment  to  a  piece  of  machinery  for  ship  quieting „ 
In  a  given  tank  and  at  a  specific  frequency  the  a  would  be  a  constant 
and  the  measurements  taken  before  and  after  treatment  would  contain 
the  same  error  and  a  true  measure  of  the  effectiveness  of  the  treatment 
would  be  given. 

A  few  comments  should  be  made  as  to  desirable  characteristics 
of  the  equipment  used.   First  an  omnidirectional  hydrophone  at  the 
frequencies  under  consideration  is  highly  desirable.   This  eliminates 
the  problem  of  orienting  the  hydrophone  and  taking  into  account  its 
directivity  ratio. 

If  a  completely  omnidirectional  hydrophone  is  not  available,, 
one  that  has  unity  directivity  ratio  in  the  horizontal  plane  is 
desirable.   This  is  so  as  to  correctly  account  for  the  direct  wave 
since  in  the  proposed  uses  of  the  technique  the  direct  wave  will 
still  be  present . 

If  the  hydrophone  is  non-directional  in  the  horizontal  plane 
care  need  not  be  taken  in  its  orientation „   Otherwise  anomalous 
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results  will  occur  due  to  the  fact  the  direct  wave  will  not  be  acting 
into  the  same  microphone  sensitivity  and  will  have  a  greater  or  lesser 
effect  depending  upon  orientation. 

Second,  the  nature  of  the  signal  to  be  measured  should  be  known c 
Possibly  it  is  not  known  exactly.   In  either  case  a  good  1/3  octave 
filter  or  band  pass  filter  should  be  used  to  ensure  that  the  power 
measured  is  of  the  desired  spectrum  in  the  event  other  unwanted 
noises  are  present. 

Third,  the  logarithmic  amplifier  should  have  a  wide  band  response 
not  only  to  handle  the  frequencies  involved  in  the  measurement  but  to 
handle  the  decay  rates  involved.   Logarithmic  amplifiers,,  in  general^ 
have  a  50  or  60  decibel  range,  the  60  db  range  being  most  desirable  „ 

The  Memoscope  has  the  potential  of  measuring  a  maximum  decay 
rate  of  5x10  db/sec„   This  is  established  on  the  basis  of  a  maximum 
sweeprate  of  10  microseconds  per  division  for  10  divisions^  50  db 
vertical  displacement  and  a  k-5   degree  slope  for  maximum  accuracy., 
Resolution  of  decay  rate  to  +  lOfo   is  quite  easily  accomplished  even 
with  little  or  no  experience  in  use  of  the  reverberation  technique 
and  evaluation  of  decay  curves. 

Lastly,  a  true  RMS  meter  similar  in  type  to  the  John  Fluke  910A 
is  recommended.   This  meter  is  unique  in  that  it  reads  true  RMS 

voltage  regardless  of  waveform  and  has  provision  for  low  or  high 

19 
damping.   When  used  with  high  damping  the  problem  of  integation 

to  get  a  voltage  measurement  is  largely  eliminated. 

19 

R.  W.  Young,  A  Brief  Guide  to  Noise  Measurement  and  Analysis^, 

USNEL,  Report  609,  May  1955 

50 


A  possible  test  setup  is  shown  in  Figure  9°   The  decade  attenuator 
and  oscillator  are  useful  in  the  calibration  of  the  logarithmic  amplifier 
Memoscope  combination. 
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APPENDIX  I 
Directivity  Ratios 

Directivity  ratios  were  determined  in  the  following  manner  since 
accurate  calibration  data  was  available  or  measured  only  for  certain 
frequencies o   From  the  beam  patterns^  horizontal  and  vertical,  the 
Directivity  Ratio  was  calculated  ,   Since  the  hydrophone  closely 
approximated  a  line  hydrophone  of  length  #,  the  equation  for  the 
pressure  field  of  a  cylindrical  line  hydrophone  is  given  in  NDRC 
Summary  as  /  Wl     v\ 


where  Y  =  angle  off  axis 
Then  since  x 
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From  Schelkunoff  the  result  of  the  integration  gives 


^  w  -  L  ^    sl  *<* 


2-* 

where  x  =  k^/2  =  Va 
This  is  then  evaluated  for  various  values  of  1/^  and  is  presented 
as  Figure  I  »  1  . 

From  the  directivity  ratios  calculated  the  curve  yielded  values 
of  J?/^  o     This  was  then  used  to  determine  an  average  lengthy  |u  for 

The  Design  and  Construction  of  Magnetostriction  Transducers, 
Summary  Technical  Report  of  Division  6  NDRC  Vol  13,  pp0  113= llU, 
Washington  DC,  19^6* 

S.  A0  Schelkunoff,  Applied  Mathematics  for  Engineers  and 
Scientists,  D.  Van  Nostrand,  19^8. 
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the  given  hydrophone ,   Then  by  use  of  this  £"  various  \/^    values 
were  calculated  for  frequencies  for  which  beam  patterns  were  not 
available.   Entering  Figure  1-1  with  J  /^   yields  the  directivity 
ratio  for  the  vertical  plane  Then  by  use  of  the  product  theorem"*, 
the  horizontal  and  vertical  directivity  ratios  were  multiplied 
resulting  in  the  directivity  ratios  given  in  Table  1=1 „ 

Table  I- 1 
Directivity  Ratios 

Frequency  Serial  Serial 

Kc  207  204 

10  1.07  1.05 

15  1.20  1.12 

20  1.355  1.25 

25  1.49  1.41 

30  1.63  1.60 

35  1.81  1.60 

40  1.99  1.77 

45  2.30  2.22 

50  2.52  2.52 


\ 


NDRC  Summary,,  op.  cit. 
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APPENDIX  II 
Hydrophone  Calibration 

Table  II- 1 

Open  circuit  receiving  response  of  LC-32^  Serial  204^,  used  as 

receiving  hydrophone  obtained  by  reciprocity  calibration  using 

pulse  technique . 

Frequency  db  re  1  volt/ubar 

Kc 

10  -102.47 

15  -103.25 

20  -104.00 

25  -103.7 

30  -105.08 

35  -105.74 

40  -106.3 

45  -107.47 

50  -107.63 
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Table  II-2 

R   Equivalent  parallel  resistance  of  LC-32  serial  207 

Frequency  Pulse  Reverberations* 

Ke  Kilohms  Kilohms 

10  1*3 ,k  35.91 

15  25-36  20.36 

20  15.47  12  Al 

25               9o75  7.81 

30               6.55  5.76 

35               ^067  ^ol2 

*+0               3A7  3c08 

^5               2.52  2o3^ 

50               1.79  lo78 

*R   for  reverberation  techniques  was  obtained  by  use  of  a  sine 

wave  oscillator  with  values  taken  at  center  frequency^,  +  2.5$^  and 

+  5$>;>  then  averaged.   For  example^,  at  lOKc  measurements  were  made 

at  10Kc,  9.5Kc,  9.75KC,  10.25Kc5  and  10.5Kc.   It  should  also  be 

commented  that  under  the  various  tank  conditions^  i.e.  no  lining^ 

six  blocks,,  12  blocks,  the  observed  R  was  within  +  3$  which  most 

likely  is  due  to  errors  in  measurement  rather  than  due  to  change 

in  R  .   It  is  expected  that  R  would  vary  under  widely  different 

tank  conditions,  but  in  the  case  investigated  the  variation  is  not 

considered  significant.   Hence,  R   is  assumed  constant  for  reverber- 

P 

ation  technique  and  is  a  source  of  a  small  error. 
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TRANSDUCER  COORDINATE 
DIAGRAM 
o 


The  spherical  coordinate  system  is  used  to  define  the  angles  in 
the  directivity  patterns  9   and  <£  shown  in  the  above  diagram  give 
the  directions  in  which  the  response  is  measured „   The  transducer 
is  placed  in  the  frame  of  reference  with  its  axis  of  symmetry 
coincident  with  6  =   0  }    its  fiducial  mark  in  the  $  =  0  plane^, 
and  its  center  at  r  =  0. 
The  two  patterns  most  frequently  measured  are; 

(1)  those  made  by  holding  <£>   constant  at  some  angle  and  rotating 
9   through  360° . 

($>  =  a  $  rotate  8) 

(2)  those  made  by  holding  9  constant  at  some  angle  and  rotating 

<£  through  360  . 

(0  =  b°|  rotate  d» 
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Figure  II-l   Horizontal  Directivity  Pattern,  Ser.  2^k,    30Kc 
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Figure  II-3   Horizontal  Directivity  Pattern,  Ser.  2Qk,    50Kc 
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Figure  II-U   Vertical  Directivity  Pattern,  Ser.  20*+,  30Kc 
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Figure  II-5   Vertical  Directivity  Pattern,  Ser.  20U,  1+OKc 
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Figure  II-6   Vertical  Directivity  Pattern,  Ser.  20U,  50Kc 
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Figure  II-7   Horizontal  Directivity  Pattern,  Ser .  207,  50Kc 
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Figure  II-8   Vertical  Directivity  Pattern,  Ser.  207,  lOKc 

1*8 


220 
140 


230" 
130 


Figure  II-9   Vertical  Directivity  Pattern,  Ser.  207,  20Kc 
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Figure  11-10   Vertical  Directivity  Pattern,  Ser .  207,  30Kc 
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Figure  11-11    Vertical  Directivity  Pattern,  Ser .  207,  ^OKc 
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